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Abstract: Based on the weak coupling expansion of gravity, we calculate the gravitational 
contributions to yukawa coupling, scalar quartic coupling as well as gauge couplings with 
general Landau-DeWitt gauge-fixing choice and a gauge preserving (of SM gauge group) 
cut off regularization scheme. We find that the results depend on the Landau-DeWitt 
gauge-fixing parameter. Based on the two loop RGE of SM couplings with one loop full 
gravitational contributions in harmonic gauge, we study the constraints on the higgs and 
top quark mass from the requirement of existing the other degenerate vacua at the Planck- 
dominated region. Our numerical calculations show that nature will not develop the other 
degenerate vacua at the Planck-dominated region with current higgs and top quark masses. 
On the other hand, requiring the existence of the other degenerate vacua at the Planck- 
dominated region will constrain the higgs and top mass to lie at approximately 130 and 
174 GeV, respectively. 
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1. Introduction 


The discovery of a 125 GeV Standard Model-like higgs boson by both the ATLAS and 
CMS collaboration SI of the Large Hadron Collider (LHC) fills the last missing piece of 
the Standard Model (SM) of particle physics. Experimental data on such scalar particle 
agree quite well with the SM predictions and no signs of new physics beyond the SM have 
been observed so far. 

Although the SM seems very successful in describing the real world, there are still 
many theoretical and aesthetical problems in SM, such as the dark matter puzzle and the 
hierarchy problem that related to the existence of a fundamental higgs scalar. It also seems 
problematic to extent the validity range of SM to Planck scale because the renormalization 
group equation (RGE) running of quartic coupling A with the current higgs and top quark 
mass will drive A negative at large field value which could lead to another local minimum at 
large field value. If such new minimum lies below our electroweak(EW) vacuum, quantum 
tunneling effects from the (false) EW vacuum to the (true) deeper vacuum could result in 
the false vacuum decay and thus the EW vacuum instability. In fact, absolute stability of 
the higgs potential is excluded at 98% C.L. for < 126GeV||3|, |, ||]. On the other hand, 
typical calculations on the tunneling rate with the central value of higgs mass indicates 
that the false EW vacuum is a metastable vacuum with a lifetime longer the age of the 
universe ip. 

A mysterious in the extrapolation of A is its slow running at high energy which is due 
to a combination of two factors: the decreasing behavior of all SM couplings at high energy 
and the nearly vanishing of at a scale of about 10^’^ to 10^® GeV. Previous SM based 
calculations indicate that the quartic coupling A and its beta function Px nearly vanish at 
the Planck scale. This may indicate the ’’multiple point criticality principle” (MPCP)0 
which assumes the other degenerate vacua at the Planck scale and was used to predict the 
top mass being 173 ± 5 GeV and a Higgs mass 135 ± 9 GeV in 1990s. However, the inputs 
of 1^] are out of date and the calculations neglect the possible gravitational contributions 
which could be important near the Planck scale. On the other hand, asymptotic safety of 
gravity]^ indicate that the MPCP may arise at a typical ’’transition scale” ktr that near 
Mpi. So it is meaningful to require the MPGP to be hold at some energy scale near the 
Mpi instead of exactly at Mpp We calculate the gravitational contributions to the beta 
functions of the standard model yukawa, gauge and quatic couplings to see the status of 
MPGP after the higgs discovery. 

The quantum effects of gravity, which is non-renormalizable by perturbative methods, 
can be studied in an effective theory approach and may play an important role near the 
Planck scale. It is interesting to note that the higgs mass had already been predicted to be 
126 GeV in the fundamental theory composed of SM and the asymptotic safety of gravity 
before LHG discovery]^]. Although gravitational effects decouple in most of the discussions 
related to standard model, such effects can change the RGE running behavior of quartic 
coupling near the Planck scale. An interesting consequence of gravitational effects is the 
asymptotic free behavior of all gange conplings near Planck scale when new power-law 
running gravitational contribntions become dominant|^. 
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Calculation of gravitational contributions to gauge couplings by |p with background 
field method was found by 0 0 to be gauge dependent and claimed that the true 
contribution vanishes. Further studies by® H, [l^ , [TsI , confirms the non-zero gravita¬ 
tional contributions to the running of gauge couplings in Q . Relevant calculations are also 
given with VilkoviskyCDeWitt effective action approach or other methods with specified 
gauge fixing condition [p!^, p0| , Our calculations are based on the traditional 

Feynman diagram methods and use a gauge invariance preserving cut-off regularization 
scheme [p3|, We should note that the weak coupling expansion of gravity used in this 
paper is not sufficient. Contributions from the non-renormalizable aspects of gravity could 
be important. The inclusion of certain higher-dimensional operators in higgs sector could 

Typical con- 


possibly modify the behavior of the potential near the Planck scale |25| , [_ 

straints from MPCP at the Planck scale on higgs and top quark masses taking into account 
the gravitational contributions was given in However, the adopted gravitational con¬ 
tributions to yukawa couplings which is very crucial in determining the UV behavior of A 
do not agree in sign with our calculations. Besides, it is more preferable to require MPCP 
in the Planck scale dominated region instead of exactly at the Planck scale. So it is very 
interesting to study the status of MPCP at the Planck scale dominated region with our 
fully independent calculations. 

This paper is organized as follows. In sec-2, we perform the calculation of gravitational 
contributions to the beta function of SM yukawa, gauge and quatic coupling with the gen¬ 
eral Landau-DeWitt gauge-fixing choice. In sec-3, we discuss the constraints from MPCP 
in the Planck scale dominated region based on our calculated gravitational contributions. 
Sec-4 contains our conclusion. 


2. Gravitional Corrections to Yukawa Coupling 


Quantum gravity is well known to be non-renormalizable. However, the quantum effects of 
gravity can be taken into account in an effective theory approach |^. Physical predictions 
for gravitational effects are justified if we interest in physics at the energy scale E < Mpi. 
Such predictions coincide with the results given by the underlying fundamental theory 
whatever its nature. 

The action S containing gauge, scalar and fermion fields can be written as 


S 


/ 


(E^XyJ—g 


k-^R- ^ 


2C 




+ E(f, + Cg + 


dg(t)du(t) - , 


H— 


( 2 . 1 ) 


where 'R' is the Ricci scalar curvature. The vierbein e“ in ( |2.1D represents the square root 
of the metric and e“^ = e°^pgP^. We can rewrite y/—g = |e^"| = e and 
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7^ = e^a'y°‘- The spin connection can be expressed in terms of vierbein as 


IJ _ (J] 


uj,, = -e 
2 




J] + 

lj,,u ' ^ ^u,pK J ) 


= ~ ~ ~ + ^e'^^e'^'^idaeuK - d^eaK)e^ .(2.2) 


The reduced gravitational coupling k = y/lGirGN ~ (1.69 x lO^^GeV) ^ are deter- 

_ 1 

mined by the Planck scale given by Mpi = ~ 1.220937 x 10^® GeV. The weak-coupling 

expansion for the Einstein gravity on Minkowski metric are given by 


dpi' — Vpi' T nh^i/ , 

gPiy = + 0{k^), 

2 

V~9 = 1 + TT^ + ~ T O(k^), (2.3) 

Z o 

where r]^u = = (1,-1,—1,-1) and h = . The vierbein can also be 

expanded in term of k 




B ^ 1 1 o 

(^pa — VPa — Vpa T '^'T'pa ~^^pp^ a T ' ' 


K 3 k-^ 

gM« = _ 1f,P^ + —h^^Php^ + 

2 8 

K 3k^ 

e^„ = eP^r]0a = S^r, - -h^r, + —h^^hpa + • 


(2.4) 


From the action, we can derive the free graviton propagator, 


pu,(7p(jp) — 


2(p2 + ie) 


Vpo-Vi'p + flppVua ^ _ 2^B^9(Tp 

1-C 




{PpPaVup + PpPpP ua H“ PuPaVfJ^p PvPpVpo 


(2.5) 


where d = 4 is the space-time dimension. In order to calculate the gravitational contribu¬ 
tions to the beta function of the yukawa coupling, we must calculate the fermion, scalar 
self energy corrections as well as the yukawa vertex corrections. 


2.1 Fermion Self Energy and The Ward-Takahashi Identity in QED 

Relevant diagrams that contribute to fermion self energy are given in fig.|^. The relevant 
lengthy Feynman diagrams are derived and given in the appendix. The contribution from 
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Figure 1: The gravitational contributions to fermion self energy. 


panel(a) is 


1 


16 Jk kP‘{k‘^ + ie) 


— [{3d -2){C-l){k-p) k‘^[{3d^ - lid + 7) + {id - 3)C] 


I 


1 


= IK 


16 Jf. k'^ + ie 

2 2 + C 

i>h, 


3d-2 


(C - 1) - {3d^ - lid + 7) - {id - 3)C 


( 2 . 6 ) 


The integral I 2 is defined as 


l2 = -i d^k 


k'^ + i 


le 


1\. 


Ibvr^ ’ 


(2.7) 


The contribution from panel(b) is 




= jD, mmp+m-k)\ ^ 


-D'^''[4}{p + k)'4){p)h{k)]V^y^„p{k) , 


f -4 [3(C + l)k^ /: + (17 - 3C)k^ p + (18C -2){k- p) 


128 
.K 


i 


[{p + ky — m? + ie] {k'^ + ie) 


= -*^(9C - 1 + 17 - 3C - 3C - 3) ^h, 

■ 2/13 3 

= + TlC) 


M--., (2.8) 

.32 32’^ - V ; 

So the sum of gravitational contributions to the self energy for fermions are given by 


iT.{^) = iK^(^ + ^C) i'^2- 
Then the renormalization constant is given as 

<^2 = ~ 1 ~ 

In our renormalizaiton, we use the known consistency!^, 

1 f d^k 

2 


k2,^,u — 


d'^k k^ki. 


(27r)4 (fe2 + A)2 


9^iu _ 1 j 

(27r)4 (A:2 + A) “ 2^^'" ' 


(2.9) 

( 2 . 10 ) 

cut off regularization rule 

( 2 . 11 ) 
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Note that the coefficient is 1/2 for quadratic divergence while be 1/4 for logarithm diver¬ 
gence. With this rule, the gauge invariance can be checked to be preserved. It is well known 
that gauge invariance requires that the fermion self energy counter term is equal to the 
fermion-antifermion-photon vertex counter term in QED. So for theoretical consistency, it 
is also important to check that the previously used renormalization procedure will preserve 
the Ward-Takahashi identity when gravity is taken into account. 

Diagrams that contribute to the gauge vertex are listed in the left panel of hgj^ which 
are calculated to be 


X X X X 







(5) 


( 6 ) 


(7) 


Figure 2: The gravitational contributions to the vertex corrections for gauge (left panel) and 
yukawa couplings (right panel). Diagrams in the second line give no contributions because they are 
logarithm divergent. 


Diagram (1): 


i{^+ /k' -I- m) 


= J Da'" + k')hip)] ,, (-»e7^) 


i{i)+ /i + m) 


K^e 


{p + k'Y — + it (p + kY — rn? -I- it 

D'Y'[Y{P + k)Y{k)h{p)\Vf,^^p„{p) , 

- 2(2C - 6)pa + (6C - 6) jhx i> 


128 

64 


/(- 2 ) 

J T) 




(p^ -I- it)[{p + kY — -I- it][(p k'Y — -I- it] ’ 

_(6 - 2C)p^7a + 8Cpa i> _ 

(p2 -I- it) [{p + k)"^ — m? + it] [{p + k')'^ — im? + it] 


= -*^(3 + 0/., 


( 2 . 12 ) 


• Diagram (2): 

-iTi = ly , 

K^el /■ 20(C - l)p^ i) - (26C-t 22 )pV 

“ 8(p2 + *e) ’ 

= (2.13) 
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Diagram (3): 


= J ^ ^TI -g -^riV^b + k)i;{k)h{p)]Vf,^^p^{k) , 

_ /■ 2[(4 - 8C)p'^ i) + (2C - 10)p^7^] 

64 Vp p^{p^ + ie)[(j) + kY - m? + ie] 

^ A f [(2 - 40 + (2C - 10)]p" 

32 "*^ Jp p'^{p'^ + ie)[{p + kY — m'^ + ie] 

= *^(C + 4)/2, (2.14) 

• Diagram (4): 

-iVY^ = f D^'^lYip + k)Y{,k)h{p)] A + m) ^ 

J p \^p ~r K) — TTl -r IC 

_ /■- 2 [( 8 C - 4)p^ ^ + (10 - 20 ^ 27 ^] 

64 Jp p2(p^ + ie)[(p + fc)2 — TO^ + le] ’ 

^ A [ [(2 - 40 + (2C - 10)]P^ 

32 "*" Jp p'^{p'^ + ie)[{p + kY - m'^ + ie] 

= *^^(C + 4)/2, (2.15) 

• Diagram (5),(6),(7) give no contributions to gauge vertex corrections because they 
are logarithm divergent. 


So the sum of the gauge vertex corrections are 

-ir = -z^r* , 

i 

= -ie-f^ — {13C+ 19)12 . (2.16) 

Then the gauge vertex counter term 5i is given by 

<^i = -^(13C + 19)/2. (2.17) 

Therefore, it is obvious from our calculation that the Ward-Takahashi identity <5i = <52 is 
satisfied with the previous choice of regularization procedure. This check is fairly non¬ 
trivial and it is justify the use of the new regularizaiton scheme. We also check that the 
Ward-Takahashi identity will be spoiled if we use the DR replacement —>■ k’^t]^'^/A for 

qudratic divergence terms. 

2.2 Yukawa Vertex Corrections and Scalar Self Energy Corrections 

We still need the gravitational corrections to yukawa vertex. Relevant diagrams that 
contribute to the yukawa vertex are listed in the right panel of fig.|^. The results are 
calculated to be: 
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Diagram (1) 


-jr(i) = / D'^'-'[^{p2)'il}{p2 + k)h{-k)\ , - —{-iY) , ^ 

Jfe P2+A — m + ie pi+/k — m + ie 

^r[^(Pi + k)'ii){pi)h{k)\D^i,^„p{k) , 

k^Y f 4C(fc^)^ — /k /k /k /k + 3(C — l)k^ Jz jk 

64 Jk [{Pi + kk — rrk + ie][{p2 + kY — mi? + ie\{k‘^ + ie) 

_ ^ r _ i2(c- i)(fc^)^ _ 

128 Jk [{pi + kY — rn? + ie][(p 2 + kY — rr? Y ieYk"^ + it) 

= *«"D^(C-1)/2 , (2.18) 

• Diagram (2) 

-*r(2) = J , 

= pcs+ 20 

8 Jk k? Y it ' 

= -iK^Y^^h , (2.19) 

• Diagram (3) 

-*r(3) = f /. _L r. ^ - -—k^n^iPi + >^)'^iPi)H-k)]'Dp^^ap{k) , 

Jk 2 piY p-mYit 

_ ^ r {2<Y -&dYk) /k /k 

16 Jk [{pi Y fc)2 — + ie](A:^ + it) ’ 

= tln^Yh , (2.20) 

O 


• Diagram (4) 

-1+6= f D^''[Y{p2)Y{P2 Y k)h{-k)] , -—-(fc) ^ 

Jk P2+ Jz-mYit 2 

k^Y r {2(Y — 6(i + 4) J: Jz 
16 Jk [iP 2 + k)"^ — rr? Y ze](fc^ + it) ’ 

= i^^kYh , ( 2 . 21 ) 

• Diagram (5) (6) (7) gives null contribution to the vertex corrections because of their 
logarithm divergence. 

So the sum of vertex correction for yukawa couplings 

7 

^Ttotal — i ^ ^ Tj , 
i=l 

= -^(27 + 290/2. (2.22) 

Then we get the renormalization constant 



The scalar self-energy corrections with the diagrams in fig. ( ^.21) had already been 
calculated in our previous work[p^]. However, the regularization procedure used in that 
calculation is the ordinary gauge non-preserving cut off scheme. So we need to reformulate 
the previous results. 




P p + k P 
b 


Figure 3: The gravitational contributions to scalar self energy. 


The scalar self-energy renormalization are calculated to be 




- / 
16 1 


2^2 


f 8d{C — l){p ■ qY — {AdC, — 8C -h — I8d + 16)p^q 


(p2 -p ieY 


2 f 2 MC - 1) - (4dC - 8C + - 18d + 16) 

7 / WTW 


(2.24) 


and 


iU{q‘^)[b] = / Ct^’'[4>{q)4>i-p - q)h]C^^[(p{p + q)(pi-q)h]Vpi,^„p{p) 

f K? A(^q^p^ 
jp 8 + if: 

• 2 2 C T 
= IK q -h . 


(2.25) 


which gives the counter term for scalar self-energy 


^0 = - 1 = ——h , 


(2.26) 


This result is different to our previous result [^| because we use the gauge preserving 
cut-off regularization scheme. 

With all the corrections at hand, we can obtain the RGE for yukawa couplings Y = 
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YoZ^z]I^Z-^ 


/3 y = ^^^Y{Yo, 

a/i /i 


= -Y 


A d 
^dA 


2 87r2 


= -A^Y{Yo,-] 

dA fx 
V + 2^<l> ~ dy 


(2.27) 


The final result is gauge dependent which is also anticipated in |31]. 

If we choose the harmonic gauge fixing (' = 1, we can obtain the gravitational contri¬ 
butions to yukawa beta function A/3y = 1. In standard model, the yukawa coupling has 
the form 


-^hfRh + h.c. , (2.28) 

so we need the replacement Y —)• yijV^ for SM yukawa couplings. The beta function from 
gravitational contribution will not be changed with such replacement. 

2.3 Gravitational Corrections To Cange And Scalar Conplings 

In order to obtain the RGE for SM, we must also include the gravitational contributions 
to gauge and scalar couplings. 




k p + k k 
b 


Figure 4: The gravitational contributions to photon self energy. 


The gravitational contributions to gauge coupling are given by the photon self-energy 
diagrams in fig. |2.3| and are calculated to be 




= -3k‘ 


2 / [(C - l)k.p{k^py + k^p^) + k^kuP^ - (C - l)(fc • p) 9 - k^iC - '^)p^iPu - P k^9 tiu] 


i 


2 p 2 _j_ 


- 3K^C{k^9f.u k^k,) I ’ 


3zk^ 


C{k^9xiu - k^iku)l2 


(2.29) 
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and 


= I + k)]G^f’^^[Ax{k)AAp + fc)]P^.,p.(p) ^ , 


= K 


I 


2Ckakxp'^ - k‘^[2Cp^ga\ - (C - ^)PaPx] - {k-p){C - ^){kxPa +Pakx) + (C - ^){k-pf‘gxoL 


— K^[(3 + C)kakx ~ (3 + C)k‘^gax] 

■ 2 3 + C , 


2[{p + ky + ze][p^ + ie] 

__ 

2[{p + ky + ie][p2 + ie] 


= IK 


-{kakx - k^gax)h , 


(2.30) 


So we can obtain 




(2.31) 


For harmonic gauge fixing ^ = 1, we can see that 


o2C-3 n 

5a = Za — 1 = k —-—12 = —K I 2 , 


(2.32) 


Therefor the gauge coupling RGE is given by 




C 2 A 2 
K A 


327r2 


(2.33) 


Our result agree with the result given in and also agrees with the traditional BFM 

(background field method) in the harmonic gauge with LORE regularization scheme |^0l| . 
Note that if we use k^k'^ —)• g^^uk '^for quadratic divergence terms, vanishing result will be 
obtained which was found in However, such replacement of DR regularization scheme 
will spoil the Ward-Takahashi identity related to gauge invariance. 

We also need the beta function for scalar quartic coupling with the new regularization 
rule k^k'^ —)• g^^^I2 for quadratic divergence terms. The relevant Eeynman diagrams can 
be seen in our previous paper [^. The non-vanishing vertex correction is given by 


2 4 


V 


{ii',ap 


(p) 


2,3 + 2C.^ 
= —K A- 1 I 2 


1 


-|- ie 


(2.34) 


The counter term 5x is given by 


5a 


-K^X 


3 + 2C 


h 


(2.35) 
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So we can obtain the beta function for quartic coupling 




= A 


d 1 + 2C Aqk^A^ 


dk 


2 167r2 

k-2A2 

= (l + 2C)Ao 


167r2 


In harmonic gauge with ^ = 1, we can obtain 

o - dX 3 
Pa = 


cZlnE' 167r2 


Xk'^E^ , 


(2.36) 

(2.37) 


We should note that the results depend on the gravitational gauge fixing choice. In 
certain circumstances, the Vilkovisky-DeWitt effective action approach which is intrinsi¬ 
cally gauge-preserving can agree with the traditional approach in harmonic gauge |Q. So 
it is preferable to fix the physical result in harmonic gauge with C = 1 such results are 
used in our subsequent calculations. 


3. MPCP Constraints With Gravitational Contributions 


Constraints from MPCP can be fairly predictive and they may reveal the existence of 
asymptotic safety of gravity. Knowing the RGE running of SM couplings, the requirement 
of MPCP near the Planck scale can be studied. On the other hand, the presence of new 
terms in the SM beta functions from gravitational effects can have important consequences. 
Such term can be dominant near the Planck mass scale and significantly change the running 
behavior of SM couplings in the UV region. 

In order to study the RGE running of quartic coupling A below the Planck scale, we 
adopt the full two-loop Standard Model beta functions|^| for A (three loop results can be 
seen in|^4|, ^), the top-yukawa couplings yt and gauge couplings gi{i = 1,2,3) in addition 
to the one-loop power-law-running gravitational contributions. Por the weak scale input, 
the following boundary conditions for RGE running are used 


as{Mz) = 0.1185 ±0.0006 
a~^{Mz) = 127.906 ±0.019 
siii^ ew{Mz) = 0.2312 ±0.0002, 
Mhiggs = 125.9 ± 0.4GeV. 
y^{Mz) = 0.0162834, 

yr{Mt) = 0 . 0102 . 


(3.1) 


which give 

a 2 {Mz) = aem(A7z)/ s\v? 6w = (29.5718)“^, 

ai{Mz) = aem{Mz)lcos^ew = (98.3341)-i. (3.2) 
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The two loop RGE running for gauge couplings in the SM are given by 


d\nE 

with 

bki — 


d bi n 

9i ~ 1 9i + 


9i 


IGvr^ 

/ 199 

9 

TT 

10 

27 

35 

10 

6 


(IhTr^)^ 


bkigl - Tr [C^fIFu + C^F^Fd + C^fIFl 

L k 


(3.3) 


44 


9 
2 

12 -26 


17 3 


10’ 2’ 


> Clfc ,cf =(^,)^,2) (3.4) 


1 3 


2 2 


3 1 


2 2 ' 


The standard GUT normalization gf = |g(y is used in the previous expressions. We 
neglect the yukawa couplings for the first two generations and keep the third generation 
contributions. The RGEs for gi are given by 


dgi ^ gi 

cZlnE IGvr^ V 10 

and 

dg2 


41 2 


1^2 


+ 


91 


(IGvr 


2'|2 


199 


27 


44 


17 


1 


dhiF IGtt^ y 6 
as well as 

dg^i 53 


52 A9 2 , , , 

' —52 H I + 


92 


din F 


IGvr^ 


I’d! + 


k^F^ 


(IGvr 


+ 


2'i2 


50 10^^ 5 10^‘ 2^^ 2^^ 


9 2|35 2 .|„2 ^2 3 2 I 2 

Y^5i + ~q92 + 1253 “ 2^* ~ ~ 2^’’ 


53 


^5? + ^92 - 265 ! - 25 ? - 25fe 


(IGtt^)^ 

The two-loop RGE for top-yukawa coupling with one loop gravitational contribution 
is given by 


d 


dlnF 


5* = 


-7 


1 


IGvr^ 

1 


^Vt + ^vlvt + vlvt - 8glyt - ^gjyt - ^9iyt + K^E^yt 


(IGtt- 

393 


2'i2 


9 


9 


-i25t^ - ^yhb - \yb + ^yhl - ^yWr - ^yt + 6A^ - i2Ayi2 - 4A5^ 


225 


+ (“oTTdi + ~rF92 + + { ^5i + tx 52 + 453 ) dfo + { ~f9i + ~f92 ) 5 


80 


16 


7 


99 


80 


16^ 


15 


15 


1187 4 9 2 2 19 2 9 23 4 22 4 

+ ■^£'1 “ + -^ 9 lgl - -^ 9 t + 95^5! - IO85I 


(3.5) 


The yukawa couplings of bottom and tau gives very small contributions to RGE of A, so 
we use here only the one-loop results 

d h 1 


din E 
d 


y = 
5 " = 


IGtt^ 

1 


^yf + ^yhb + 5rdb - - ^525?) - 853^6 


^3iq2 iq2 92 92 

2dr + ^Vt Vt + ^ybyr — ^5l — ^52 


dlnE’" IGtt^ 

The normalization of higgs potential in our study is given as ^ 


m 


A 


V = -+ -h^ . 


(3.6) 


(3.7) 


^ There is a factor 1/2 difference for A with respect to the convention used in 
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The improved two-loop RGE for A with one loop gravitational contribution is given by 


tA = 




+ 


(IGvr 


■ 2\2 


Pi, 


(3.8) 


dlnE' 

within which the one loop part is 

Px = 24A^ - 2 (tOO^^ 

~ 2(3?/^ + ^Vb + vt) + SAk^E^ , (3.9) 

and the two loop part[ 


IS 


^9 .108 n n, . /73 4 117 o o 

Px — —312A -I- i~^9i + 1085'2)A — ^9i92 


1887 ^4 
200 ‘ 


9t]X 


+ 2 


1 /305 g 289 


■92 


40 


4 2 
-9291 


1677 


-9291 


1491 

1000 


9'i 


- + 2/fe) - {‘^yf -yt + ^yt) - + ^yj + y 


+ lOA 


2 


/^11'2|0 2j_o 2\ 2|/^l2|0 2|o2\ 2|/^3 2 3 2\ 2 
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The inputs for RGE, namely the MS gauge coupling, higgs self coupling and top yukawa 
coupling at the pole top mass, can be obtained with Q 


X{mt) 

yt{mt) 


gsim) 


0.12577 -h 0.00205 (_ 125 ) _ 0.00004 (- 

0.93587 -h 0.00557 (- 173.15) - 0.00003 ( 

VGeV / V GeV 

'as(Mz) - 0.1184\ 

^ ^ -j +0.00200*;,, 

' ' rnt 


0.00041 


1.1645 + 0.0031 


0.0007 

as(Mz) -0.1184\ 


0.0007 




- 0.00046 


GeV 


173 . 15 ) ±0.00140*^, 

- 125 ) 


- 173 . 15 ) . (3.11) 


An illustration of full gravitational contributions to the RGE running of A are shown in 
fig.5. We can see that the gravitational contributions greatly modify the RGE trajectory 
near the Planck scale. 

We scan the parameter space of {Mh, Mt) to check the points which can lead to MPGP 
at some transition scale near Mpi. We dehne the Planck-scale dominated region to be: 
8.0 X 10^’^GeV < E < Mpp We use the following input: 


• The existence of degenerated vacua near the Planck scale can be satisfied if |/3 a(E)| < 
1.0 X 10“^, |A(E)| < 1.0 X lO”"^ and also d?\{E)/dE'^ > 0. 

Our scan indicate that the MPGP condition at the Planck scale dominated region 
can not be satisfied with current top quark mass m* G (173.21 ± 1.22) GeV or higgs 
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Figure 5: The left panel denotes the gravitational contributions to the RGE running of various 
couplings: A, j/t,5i, 52 , ffa- The right panel compares the RGE of quartic coupling A near the Planck 
scale with and without gravitational couplings. 


mass ruh G (125.09±0.32) GeV. The numerical results for the MPCP constraints with 
current measured top or higgs values are shown in fig.|^ We can see that MPCP could 
be satisfied at approximately 10^^ —10^^'® GeV for the current measured values of top 
and higgs masses. Gravitational effects are negligible in finding the large-field-value 
degenerate vacua located in this region. 
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Figure 6 : The points with current measured higgs and top quark masses that satisfy the MPCP 
conditions with the setting \P\{E)\ < 1.0 x 10“^ and |A(i?)| < 1.0 x 10“'^. 

• We scan the parameter space for higgs range nih G [120,135] GeV and top range 
mt G [165,180] GeV. We keep those points that correspond to new degenerate vacua 
in the range 8.0 x 10^^ < E < Mpi (the Planck-scale dominated region). The results 
can be seen in the left panel of fig.|^ From the the left panel, we can see that the higgs 
mass is constrained to lie in the range [129.0,130.2] GeV while the top quark mass 
in the range [173.8,174.4] GeV. Our scan also indicates that the highest degenerated 
vacua lies at approximately 2.2 x 10^® GeV. However, if we relax the MPCP criterion 
to |A(i?)j < l.Ox 10“^, it is possible to obtain the other degenerate vacua at the Planck 
scale with multiple choices of top quark and higgs mass, for example, mt = 175.10 
and mfi = 132.14. Constraints for MPCP in the Planck-scale dominated region also 
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give an upper bounds on higgs and top quark masses as rrih < 134.9 and mt < 176.34, 
respectively. 
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Figure 7: The results of the scan in the parameter space mtiggs & [120,135] and mtop & [165,180] 
to seek the degenerated vacua in the Planck-dominated range (8 x lO^’^ < E < Mpi). Each panel 
denotes the projection of the surviving triples {E,mt,mh) onto the corresponding plane. 


• It is instructive to compare our scan result with the case without gravitational con¬ 
tributions. Without gravitational effects, the bounds for higgs mass and top quark 
masses can be released. We can see from fig.§ that the higgs mass within [129,135] 
GeV and top quark within [173,176] GeV can lead to degenerated vacua in the Planck- 
dominated region. 
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Figure 8: Dependence on the value of as{Mz) for higgs and top quark masses that can lead to the 
other degenerated vacua in the Planck-dominated range (8 x lO^"^ < E < Mpi) without(the upper 
panel) and with (the lower panel) gravitational effects. 
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• Our scan also show that the results are sensitive to the value of as{Mz)- We show the 
results of our scan with different choice of as{Mz) in figj^. Larger as{Mz) requires 
larger mtop and ruhiggs masses. 


4. Conclusions 

Based on the weak coupling expansion of gravity, we calculate the gravitational contribu¬ 
tions to yukawa coupling, scalar quartic coupling as well as gauge couplings with general 
Landau-DeWitt gauge-fixing choice and a gauge preserving (of SM gauge group) cut off 
regularization scheme. We find that the results depend on the Landau-DeWitt gauge-fixing 
parameter. Based on the two loop RGE of SM couplings with one loop full gravitational 
contributions in harmonic gauge, we study the constraints on the higgs and top quark 
mass from the requirement of existing the other degenerate vacua at the Planck-dominated 
region. Our numerical calculations show that nature will not develop the other degenerate 
vacua at the Planck-dominated region with current higgs and top quark masses. On the 
other hand, requiring the existence of the other degenerate vacua at the Planck-dominated 
region will constrain the higgs and top mass to lie at approximately 130 and 174 GeV, 
respectively. 
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Appendix A: Feynman rules 

To simply the lengthy expressions in certain Feynman rules, we define the following ex¬ 
pressions: 


Cfj,u^pcr — '^fip'^ucr T 'Hi^cr^vp 'npiL'Vpo'' 

Cp, ,pcr(p) — P^Pp'Hva PpPaPvp + PuPpPpa + PuPaPpp 


^ 2 ) — Ppuk\(jk2p 


Ppak\ifk2p Pppk\o-k2i/ Ppak\pk2i/ (/^ ^ 


^ph>,p<j\oJT 2 “1“ PaLjC-\- PpujCPuujCp,r^pCF PujtC + (w O* t)] 


kkpvp(TUJT{kl’)k2^ — 


,pT^l(T^2a; “1“ Cpy^pijjk\q-k2(j C^p,L',a;T kipk2u + (/o ^ < 7 ) 


(/r,i/) o (p,fT) 
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IflUptJUJT {ki, k2) + 

{CaUJ,VT 'n(TT'UuUj)k\^k2p “1“ {CuLOjfTT ^(TOJ ^2/i “1“ (M 


+ (po a) 


We can derive the Feynman rules from the weak gravity expansion of metric tensor 
on Minkowski background. The vierbein and spin connections can also be expanded ac¬ 
cordingly. After lengthy algebraic manipulations, we can obtain the relevant Feynman 
rules: 




+ ^2 - 2m) + [(Cpe,,„a - \vvpVcra - {kl + fc2)p 

-^inupla + 'Ji/o-7p)(Pl - P2)p + \{l^l - j/2)(’7pCT'^pp - tlup'lalp) + ((i ■<-»■ 1^)] 
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